A microscopic imaging method of phase contrast and differential interference contrast based on the transport of intensity equation

ABSTRACT

The invention claims a microscopic imaging method of phase contrast (PC) and differential interference contrast (DIC) based on the transport of intensity equation (TIE). Firstly, three intensity images are captured along the optical axis; secondly, TIE is solved by deconvolution to obtain the quantitative phase; then, the intensity image under the DIC imaging mode is obtained according to the DIC imaging principle; finally, the corresponding phase image of PC imaging mode is obtained according to the PC imaging principle. The proposed approach can endow the bright-field microscope with the ability to realize PC and DIC imaging without complex modification of the traditional bright-field microscope. In other words, this method only needs to use the traditional bright-field microscope without adding any complex hardware. Through the PC and DIC algorithms, this method has the advantages of quantitative, high-speed, low-cost, simple structure, and less external interference. In addition, it has the same imaging performance as the phase contrast microscope and differential interference contrast microscope, which are expensive, complex-structure, and has strict environmental conditions.

FIELD OF THE INVENTION

The invention belongs to the optical microscopic measurement, imaging technology, especially the microscopic imaging method of phase contrast and differential interference contrast based on the transport of intensity equation.

BACKGROUND OF THE INVENTION

In the field of biomedical microscopy imaging, most samples are phase objects with a uniform distribution of amplitude transmittance and an uneven distribution of refractive index or thickness. Therefore, when light waves pass through a phase object, the wavelength and amplitude change very little, while the phase changes dramatically (Schmalz J A, Gureyev T E, Paganin D M, et al. Phase retrieval using radiation and matter-wave fields: Validity of Teague's method for solution of the transport-of-intensity equation[J]. Physical Review A, 2011, 84(2): 023808). However, the human eye and optical detectors are unable to observe this phase difference with important information. To address this challenge, academics usually take the approach of staining the sample. The different affinities of different intracellular components for different chemical or fluorescent dyes are exploited to create a sufficiently large contrast in intensity or to generate different spectra for the purpose of cell imaging. Fluorescence microscopy and laser confocal microscopy are widely used to display the structural and functional properties of cells by selectively labeling specific molecules within the cells. However, both traditional dye staining methods and fluorescent labeling staining methods are not conducive to cellular research in the biomedical field because the labeling means are not able to observe living cells for a long time due to the damage to the cells (Sander P, Nehorai A. Deconvolution methods for 3-D fluorescence microscopy images[J]. IEEE Signal Processing Magazine, 2006, 23(3): 32-45).

In the study of dynamic microscopic imaging of living cells, label-free imaging is an effective imaging modality. It can be divided into phase visualization methods and phase measurement methods. The most common phase visualization methods are Phase Contrast (PC) microscopy and Differential Interference Contrast (DIC) microscopy. In 1935, the Dutch scientist Zernike invented phase contrast microscopy, which transformed the optical path length through different parts of the object into a difference in amplitude (light intensity), allowing the observation of unstained specimens by obtaining phase images with high contrast (Zuo Chao, Chen Qian, Sun Jiasong. Non-interferometric phase recovery and quantitative phase microscopic imaging based on the light intensity transport equation: a review of the literature and recent advances [1]. (China Laser, 2016, 43(6): 0609002.). However, this method can obtain qualitative phase images, does not involve quantitative analysis, and requires the use of a special condenser with annular diaphragm and phase contrast objective lens with phase plate, which are complex and expensive. 2015 Laura Waller's group in the United States proposed a digital phase contrast technique based on array light sources with structured light illumination, which uses LED array to obtain microscopic images with different illumination directions. The digital phase contrast algorithm is used to synthesize digital phase contrast microscopic images. This method has a simple structure and can be used to obtain quantitative phase contrast images quickly. However, due to the limitation of the numerical aperture of the light source, it is not effective for imaging phase objects (cells, etc.) (Tian L, Waller L. Quantitative differential phase contrast imaging in an LED array microscope [J]. Optics express, 2015, 23(9): 11394-11403). 2017 Guangdong OPPMANDI Technology Co., Ltd. proposed a dual-channel structured light digital phase contrast microscope imaging system (CN206920698U), which reduces external perturbation by simultaneous acquisition of images through an image acquisition module composed of two cameras, with less interference from the outside world and can obtain the phase contrast image in real-time. However, this method requires the addition of an unpolarized spectroscopic prism and an additional camera, which makes the structure complicated and the cost increased.

DIC microscopy is a wavefront shear interference technique that uses the principle of polarized light interference. The phase change of the sample can be expressed in the form of intensity, where the areas with large phase gradients have high contrast in the field of view and show a “pseudo-stereoscopic” effect (Fu D, Oh S, Choi W, Yamauchi T, Dorn A, Yaqoob Z, et al. Quantitative DIC microscopy using an off-axis self-interference approach. Opt. Lett. 2010; 35:2370-2372). However, this method only supports qualitative analysis of phase images and requires the use of complex experimental equipment: two sets of polarizers and polarization beam splitting prisms (Wollaston prism or Nomarski prism). In 2012, Zeng Shaoqun of Huazhong University of Science and Technology and others proposed a super-resolution differential interference phase contrast microscopy imaging system and microscopy imaging method (CN102998789A), which solved the problem of low imaging contrast of structured light and scattered light imaging, thus realizing the high-contrast super-resolution imaging of biological sample without processing. However, the spatial light modulator used is expensive, and its refresh speed limits the frame rate of image acquisition. It cannot be applied to DIC imaging in the high-speed changing process. In addition, there is the problem that the polarization state of light affects the modulation effect of light by the spatial light modulator

Therefore, how to achieve quantitative, high-speed, and low-cost phase contrast and differential interference contrast imaging in a bright-field microscope without additional complex devices has become a technical challenge in the observation of label-free samples by computational microscopy imaging.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide a microscopic imaging method of phase contrast (PC) and differential interference contrast (DIC) based on the transport of intensity equation (TIE), which can realize phase contrast and differential interference contrast imaging by bright-field microscopy without complex modification of the conventional bright-field microscope.

The technical solution to achieve the purpose of the present invention is:

Compared with the prior art, the invention has significant advantages: (1) Simple structure. It can be directly realized on the traditional bright-field microscope without adding the complicated hardware required for traditional PC imaging and DIC imaging. The proposed method can achieve the same imaging effect as traditional PC imaging and DIC imaging (expensive, complex, and demanding in environmental conditions), which not only reduces the complexity of the device but also ensures the imaging quality. (2) Non-interference. It can work well under partially coherent illumination, so the working environment is not demanding and the application field is wider. (3) Low cost. Only the traditional bright-field microscope does not need to use additional complicated devices, which solves the problem of expensive components in the existing phase contrast imaging system, and avoids the problem of the spatial light modulator being sensitive to the polarization of light.

The invention is further described in detail below with reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic diagram of the optical path of the microscopic imaging microscope based on the transport of intensity equation for phase contrast and differential interference contrast.

FIG. 2 is a flowchart of the microscopic imaging method of phase contrast and differential interference contrast based on the transport of intensity equation.

FIG. 3 shows the schematic diagram of different phase contrast rings in the phase contrast method. FIG. 3(a) is the phase contrast ring under annular phase contrast; FIG. 3(b) represents the relationship between the amplitude transmittance (T) and the ring radius (r) of the phase contrast ring half-section under annular phase contrast; where FIG. 3(c) shows the relationship between the phase (P) and the ring radius (r) of the phase contrast ring half-section under annular phase contrast (or apodization phase contrast) (positive phase contrast case); where FIG. 3(d) displays the phase contrast ring under apodization annular phase contrast FIG. 3(e) shows the relationship between the amplitude transmittance (T) and the ring radius (r) of the phase difference ring half-section under apodization phase contrast; FIG. 3(f) represents the relationship between the phase (P) and the ring radius (r) of the phase difference ring half-section under annular phase contrast (or apodization phase contrast) (negative phase contrast case).

FIG. 4 displays the simulated images of the microscopic imaging method of phase contrast and differential interference contrast based on the transport of intensity equation. FIG. 4(a 1) shows the phase difference plate under the ring phase contrast; FIG. 4(a 2) represents the phase difference plate under the apodization phase contrast; FIG. 4(a 3) shows the Amplitude Transfer Function (ATF) of the ring phase contrast; FIG. 4(a 4) is the ATF of the apodization phase contrast; FIG. 4(b 1) shows the Phase Transfer Function (PTF) of the ring positive phase contrast; FIG. 4(b 2) is the PTF of the apodization positive phase contrast; FIG. 4(b 3) displays the PTF of the annular negative phase contrast; FIG. 4(b 4) shows the PTF of the apodization negative phase contrast; FIG. 4(c 1) shows the PTF of the annular phase contrast under annular illumination with the numerical aperture of the objective lens of 0.4; FIG. 4(c 2) represents the annular phase contrast PTF of the annular illumination with the numerical aperture of the objective lens of 0.65; FIG. 4(c 3) shows the annular phase contrast PTF under annular illumination with a numerical aperture of the objective lens of 0.75; FIG. 4(c 4) is the annular phase contrast PTF under annular illumination with a numerical aperture of the objective lens of 1.25; FIG. 4(d 1) is the apodization phase contrast PTF under annular illumination with a numerical aperture of the objective lens of 0.4; FIG. 4(d 2) represents the apodization phase contrast PTF under annular illumination with the objective lens numerical aperture at 0.65; FIG. 4 (d 3) shows the apodization phase contrast PTF under annular illumination at a numerical aperture of the objective lens at 0.75; FIG. 4 (d 4) expresses the apodization phase contrast PTF under annular illumination at a numerical aperture of the objective lens at 1.25.

FIG. 5 shows a schematic diagram of the HeLa cell results of the present invention. FIG. 5 (a) displays the quantitative phase image obtained from the transport of intensity equation (TIE); FIG. 5 (b) shows the annular positive phase contrast light intensity image; FIG. 5 (c) expresses the apodization positive phase contrast light intensity image; FIG. 5 (d) shows the differential interference contrast (DIC) light intensity image; FIG. 5 (e) represents the annular negative phase contrast light intensity image; FIG. 5 (f) displays the apodization negative phase contrast light intensity image.

DESCRIPTION OF THE PREFERRED EMBODIMENT

As shown in FIG. 1, the present invention is based on a microscopic imaging system equipped with an adjustable lens module with a standard 4F relay system, the actual hardware platform of which is an inverted Olympus IX71 microscope (it is composed of condenser 1, aperture stop 2, condenser lens 3, objective lens 5, and imaging tube lens 6 in FIG. 1). The microscope was equipped with an Olympus camera 10 (Olympus UC50, resolution 2588 pixels×1960 pixels, 3.4 μm/pixel) for image acquisition and an electronically controlled zoom lens (EL-C-10-30-VISLD, Optotune AG) module 11 for precise zooming. This system is controlled by software via the USB interface and ensures the synchronization of the camera 10 acquisition and the electronically controlled zoom lens module 11 zoom in the direction perpendicular to the focal plane (z-axis). The image stack is acquired by camera 10 after passing through a long working distance planar semi-complex achromatic objective 5 (Olympus, LUCPLFLN 40×, NA 0.6). It can be seen from FIG. 1 that the back focal plane of the objective lens 5 coincides with the front focal plane of the imaging tube lens 6, and the imaging plane of the camera 10 is placed at the back focal plane position of the 4F relay system behind the imaging tube 6 (composed of lenses 8, 9 in FIG. 1). During imaging, the sample 4 to be tested on the stage is adjusted to the position of the front focal plane of the objective lens 5 to form an infinity-corrected imaging system. An electronically controlled zoom lens module 11 is placed on the spectrum surface of the standard 4F relay system. The light is condensed by the condenser 1 and becomes part of the coherent light that illuminates the sample 4 to be tested. The sample 4 to be tested is placed on the stage. The light passes through the sample 4, passes through the imaging tube lens 6, and then converges through the electronically controlled zoom. The lens module 11 is modulated to illuminate the imaging plane 7 of the camera 10, and a series of intensity images can be collected.

With reference to FIG. 2, the present invention utilizes the above-mentioned device to realize the phase contrast (PC) and differential interference contrast (DIC) based on the transport of intensity equation (TIE), including the following four steps:

Step 1, three intensity images are acquired along the optical axis, and the process is as follows: three intensity images are acquired along the optical axis: positive defocus image I_(z-P), infocus image I₀, and negative defocus image I_(z-N). The first-order axial differentiation ∂I/∂z of the intensity I is obtained by using numerical difference estimation

$\frac{\partial I}{\partial z} = \frac{I_{z - P} - I_{z - N}}{dz}$

where dz is the distance between the positive defocus image I_(z-P) and the negative defocus image I_(z-N) in the Z axial direction, and the value is generally 5×10⁻⁶ m.

Step 2, obtain the quantitative phase image ϕ based on the TIE. The procedure is as follows: solve the TIE using the first-order axial differentiation ova, obtained in Step 1

${{- k}\frac{\partial I}{\partial z}} = {\nabla{\cdot \left( {I{\nabla\phi}} \right)}}$

obtain the quantitative phase image ϕ

$\phi = {{- k}{{\nabla^{2}\nabla} \cdot \left( {I^{- 1}{\nabla{\nabla^{- 2}\frac{\partial I}{\partial z}}}} \right)}}$

where k is the wave number, it can be expressed as k=2π/λ (λ is the wavelength). And ∇ is the gradient operator, · denotes the dot product, ∇⁻² is the inverse Laplacian operator.

Step 3, solve TIE in the differential interference contrast (DIC) imaging mode. The specific implementation process is:

(1) Firstly, a beam of light passes through a Wollaston prism and divides into two coherent light fields with Δx, Δy shear;

(2) Any one of them produces another phase shift of π/2, the other one remains unchanged;

(3) Finally, the two beams interfered and superimposed to form an intensity interferogram, that is, the final DIC image obtained. The specific process is as follows: First, the resulting object optical field U(x,y) is:

U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)]

where I₀ is the focus image obtained in Step 1, ϕ(x,y) is the ϕ obtained in Step 2, then the light field U_(Δ)(x,y) after shearing Δx,Δy and phase shifting π/2 is

U _(Δ)(x,y)=U(x+Δx,y+Δy)=√{square root over (I/(x+Δx,y+Δy))}exp{i[ϕ(x+Δx,y+Δy)+π/2]}

then the intensity image DIC(x,y) formed by U(x,y) and U_(Δ)(x,y) is

DIC(x,y)=|U(x,y)+U _(Δ)(x,y)|²=|√{square root over (I ₀)} exp[iϕ(x,y)]+√{square root over (I/(x+Δx,y+Δy))}exp{i[ϕ(x+Δx,y+Δy)+π/2]}|²

where Δx,Δy is integer N, its value is used to adjust the amount of shear and the positive and negative is used to adjust the direction of shear; DIC(x,y) is the obtained DIC intensity image. The solution for the special case of DIC imaging is as follows: {circle around (1)} For the case of uniform intensity there has I₀=1, then

U(x,y)=exp[iϕ(x,y)]

U _(Δ)(x,y)=U(x+Δx,y+Δy)=exp{i[ϕ(x+Δx,y+Δy)+π/2]}

DIC(x,y)=|U(x,y)+U _(Δ)(x,y)|²=|exp[iϕ(x,y)]+exp{i[ϕ(x+Δy,y+Δy)+π/2]}|²

{circle around (2)} For weakly phased objects for uneven intensity, i.e., ϕ(x,y)□ 1, so

U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)]=√{square root over (I ₀)}[1+iϕ(x,y)]

DIC(x,y)=I ₀{1+[ϕ(x,y)−ϕ(x+Δx,y+Δy)]}

{circle around (3)} For weakly phased objects for uniform intensity, thus

  U(x, y) = exp [i ϕ(x, y)] = 1 + i ϕ(x, y) U_(Δ)(x, y) = U(x + Δx, y + Δ y) = exp {i[ϕ(x + Δx, y + Δy) + π/2]} = exp (i π/2)[1 + i ϕ(x + Δx, y + Δy)] = i[1 + i ϕ(x + Δx, y + Δy)] = i − ϕ(x + Δx, y + Δy) ${DI{C\left( {x,y} \right)}} = {{{{U\left( {x,y} \right)} + {U_{\Delta}\left( {x,y} \right)}}}^{2} = {2 + {2\left\lbrack {{\phi\left( {x,y} \right)} - {\phi\left( {{x + {\Delta x}},{y + {\Delta y}}} \right)}} \right\rbrack} + {\underset{\underset{\approx 0}{︸}}{{\phi^{2}\left( {{x + {\Delta x}},{y + {\Delta y}}} \right)} + {\phi^{2}\left( {x,y} \right)}}{\bullet 1}} + \left\lbrack {{\phi\left( {x,y} \right)} - {\phi\left( {{x + {\Delta x}},{y + {\Delta y}}} \right)}} \right\rbrack}}$

In this case, ϕ²(x+Δx,y+Δy)+ϕ²(x,y) can be approximated as zero, so the DIC intensity image obtained in this case is

DIC(x,y)=1+[ϕ(x,y)−ϕ(x+Δx,y+Δy)]

{circle around (4)} For the case of slowly varying intensity √{square root over (I₀)}≈√{square root over (I(x+Δx,y+Δy))}, there is

DIC(x,y)=|√{square root over (I(x+Δx,y+Δy))}{exp[iϕ(x,y)]+exp{i[ϕ(x+Δx,y+Δy)+π/2]}}|².

(4) In practice, however, it is not only limited to integer shear quantities but also fractional shear quantities can be achieved. Then the solution, in this case, must be performed in the frequency domain:

U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)]

U _(Δ)(x,y)=exp(iπ/2)FFT⁻¹ {Û(u _(x) ,u _(y))exp[i2π(u _(x) Δx+u _(y) Δy)]}

DIC(x,y)=|U(x,y)+U _(Δ)(x,y)|²

The above method is a sub-pixel translation based on the Fourier domain, where Û(u_(x),u_(y))=FFT{U(x,y)} is the Fourier transform of, (u_(x),u_(y)) is the Fourier-domain coordinate corresponding to the space domain coordinate (x,y) (i.e., the frequency-domain coordinate). The value of the shear Δx,Δy is adjustable. When the shear Δx,Δy is integer, it can be solved directly in the space domain, and when the shear Δx,Δy is fractional number, it needs to be converted to the frequency domain. Meanwhile, the magnitude of the shear Δx,Δy determines the convexity of the object in the DIC image, which is often within the range of (1,10); and the positive or negative shear value Δy,Δy determines the shading direction of the object in the DIC image.

This step results in an intensity image in the DIC imaging mode DIC(x,y).

Step 4, solve the phase transfer function in the phase contrast (PC) imaging mode, and solve the intensity image in the PC imaging mode by the annular phase contrast method or the apodization phase contrast method.

The specific implementation process is: using the quantitative phase image obtained in Step 2, based on the phase contrast imaging theory to obtain the intensity image I_(PC-A)(u) obtained by the annular phase contrast method and the intensity image I_(PC-B)(u) obtained by the apodization phase contrast method. The first method is the annular phase contrast method: compared with the traditional phase contrast method, the difference is that the illumination is annular (instead of plane waves) and the phase contrast plate in the imaging system is annular (instead of a point in the middle) [see FIG. 3(a), FIG. 3(b), FIG. 3(c), and FIG. 3(f) for the specific forms]. The intensity and phase of the object are obtained by phase recovery and the resulting object optical field U(x,y)

U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)]

Thus,

I(u) = Bδ(u) + A(u)ATF(u) + iΦ(u)PTF(u) B = ∫∫L(u_(j))P(u_(j))²d²u_(j) ${{AT}{F(u)}} = {{\frac{- {\int{\int{{{L\left( u_{j} \right)}\left\lbrack {{{P^{*}\left( u_{j} \right)}{P\left( {u + u_{j}} \right)}} + {{P\left( u_{j} \right)}{P^{*}\left( {u - u_{j}} \right)}}} \right\rbrack}d^{2}u_{j}}}}}{\int{\int{{L\left( u_{j} \right)}{{P\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}PT{F(u)}} = \frac{\int{\int{{{L\left( u_{j} \right)}\left\lbrack {{{P^{*}\left( u_{j} \right)}{P\left( {u + u_{j}} \right)}} - {{P\left( u_{j} \right)}{P^{*}\left( {u - u_{j}} \right)}}} \right\rbrack}d^{2}u_{j}}}}{\int{\int{{L\left( u_{j} \right)}{{P\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}}$

where I₀ is the focus image obtained in Step 1, ϕ(x,y) is the phase image ϕ obtained in Step 2, u denotes the frequency coordinate in Fourier space; I(u) is the intensity spectrum distribution; B indicates background items; A(u) is the amplitude spectrum distribution of the sample; ATF(u) denotes the amplitude transfer function; Φ(u) is the phase spectrum distribution of the sample; PTF(u) denotes the phase transfer function; u_(j) is the corresponding frequency shift vector caused by a single tilted illumination, j indicates different angles of illumination; L(u_(j)) denotes the intensity distribution of the illumination source in the Fourier domain; P(u) is the pupil function of the objective lens; * denotes the conjugate operation. The above expressions are general expressions for the background term, amplitude transfer function and phase transfer function under complex illumination and are applicable to all illumination modes and pupil functions; therefore, in the case of annular phase contrast, the light source expression for the annular illumination mode L(u) is

${L(u)} = \left\{ \begin{matrix} {1,} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {0,} & {\ {others}} \end{matrix} \right.$

The pupil function under the annular phase contrast P_(PC-A)(u) is

${P_{{PC} - A}(u)} = \left\{ {{\begin{matrix} {1,} & {k_{1}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{0}^{2}} \\ {{0.25 \times {\exp\left( {i\frac{\pi}{2}} \right)}},} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {1,} & {{u_{x}^{2} + u_{y}^{2}} < k_{2}^{2}} \\ {0,} & {others} \end{matrix}k_{i}} = {N\;{A_{i}/{\lambda\left( {{i = 0},1,2} \right)}}}} \right.$

where (u_(x),u_(y)) is the Fourier-domain coordinate corresponding to the space domain coordinate (x,y), λ is the wavelength, k₀ is the radius of the corresponding objective limit, k₁ is the outer diameter of the phase contrast ring, k2 is the inner diameter of the phase contrast ring. In the set of simulation experiments shown in FIGS. 4(a 1), 4(a 3), 4(b 1), 4(b 3), and 4(c 1), the NA₁=0.28, NA₂=0.25, NA₀=0.4 are used. The ring illumination integration region corresponds to the phase difference ring (inner diameter is k₂, outer diameter is k₁). Then the annular phase contrast of ATF_(PC-A)(u) and PTF_(PC-A)(u) can be expressed as

${{AT{F_{{PC} - A}(u)}} = \frac{\begin{matrix} {- {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{PC} - A}^{*}\left( u_{j} \right)}P_{{PC} - A}\left( {u + u_{j}} \right)} +} \right.}}}} \\ {\left. {{P_{{PC} - A}\left( u_{j} \right)}{P_{{PC} - A}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{PC} - A}\left( u_{j} \right)}}^{2}d^{2}u}}}}{{PT{F_{{PC} - A}(u)}} = \frac{\begin{matrix} {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{PC} - A}^{*}\left( u_{j} \right)}P_{{PC} - A}\left( {u + u_{j}} \right)} -} \right.}}} \\ {\left. {{P_{{PC} - A}\left( u_{j} \right)}{P_{{PC} - A}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{PC} - A}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}}$

The corresponding with ATF_(PC-A)(u) and PTF_(PC-A)(u) results are shown in FIG. 4(a 3) and FIG. 4(b 1). Under complex partially coherent illumination like annular phase contrast, the intensity distribution of the acquired image can be regarded as a linear superposition of the intensity at different angles, so we decompose the complex illumination into illumination at different angles for analysis. At this point, the annular illumination can be represented by the Fourier translation technique, so the Fourier transform of the complex amplitude on the camera plane satisfies

W _(j)(u)=√{square root over (L(u _(j)))}[δ(u−u _(j))−A(u−u _(j))+iΦ(u−u _(j))]P(u)

Thus, the spectral distribution of the intensity images acquired by the camera under single-angle tilt illumination can be written as a convolution of W_(j)(u) and its complex conjugate w_(j)*(u):

(u)=W _(j)(u)⊗W _(j)*(u)=L(u _(j))δ(u)|P(u _(j))|² −L(u _(j))A(u)[P*(u _(j))P(u+u _(j))+P(u _(j))P*(u−u _(j))]+iL(u _(j))Φ(u)[P*(u _(j))P(u+u _(j))−P(u _(j))P*(u−u _(j))]

In order to linearize the phase recovery, the cross-convolution term in the above equation is neglected. When the sample is partially illuminated by spatially coherent light, the resulting annular phase contrast intensity spectrum I_(PC-A)(u) can be written as the sum of the intensity spectrum at multiple illumination angles

I _(PC-A)(u)=B _(PC-A)δ(u)+A(u)ATF_(PC-A)(u)±iΦ(u)PTF_(PC-A)(u)

For the phase difference plate, the smaller the size, the higher the contrast of the low-frequency component; for the attenuation rate, it can be used to control the ratio of the low-frequency component to the high-frequency component—usually the larger the value taken the greater the attenuation and the higher the phase contrast (90%); for the phase shift amount, there are two cases of positive and negative phase contrast, the difference lies in the phase in the phase difference plate that the delay of the former is +π/2 and the delay of the latter is −π/2. The algorithm can obtain different annular phase contrast images by changing the value of the numerical aperture of the objective NA₀, where the simulation results of changing the numerical aperture of the objective NA₀ are shown in FIG. 4(c 1) to FIG. 4(c 4).

Another method based on phase contrast imaging theory is the apodization phase contrast method: the difference compared to the annular phase contrast method is the use of an apodization phase contrast plate. By using selective amplitude filters in the phase contrast plate located in the objective lens in the back focal plane adjacent to the phase film, the effect of reducing halos and improving the contrast of the sample can be achieved. These amplitude filters consist of neutral density filter films, and the films are applied to the phase contrast plate around the phase film. The transmittance of the phase shift ring in the latest variant phase difference plate is about 25%, while the pair of adjacent rings around the phase shift ring in the variant phase difference plate has neutral density, and the transmittance varies uniformly from 25% to 100% on both sides [see FIG. 3(d), FIG. 3(e), FIG. 3(c), and FIG. 3(f) for the specific form], and the width of the phase film is the same in both plates. The solution method is similar to the circular phase contrast, the difference is that the former is a step filter and the latter is a trapezoidal filter, i.e., the pupil function P_(PC-B)(u) is different. It is expressed by the equation

${P_{{PC} - B}(u)} = \left\{ {{\begin{matrix} {1,} & {k_{1}^{\prime 2} < {u_{x}^{2} + u_{y}^{2}} < k_{0}^{2}} \\ {{{\frac{3}{4\left( {k_{2}^{\prime} - k_{1}} \right)}u_{x}} + \left\lbrack {1 - {\frac{3}{4\left( {k_{1}^{\prime} - k_{1}} \right)}k_{1}}} \right\rbrack}\ ,} & {k_{1}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{\prime 2}} \\ {{0.25 \times {\exp\left( {i\frac{\pi}{2}} \right)}}\ ,} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {{{\frac{- 3}{4\left( {k_{2} - k_{2}^{\prime}} \right)}u_{x}} + \left\lbrack {1 - {\frac{3}{4\left( {k_{2} - k_{2}^{\prime}} \right)}k_{2}}} \right\rbrack}\ ,} & {k_{2}^{\prime 2} < {u_{x}^{2} + u_{y}^{2}} < k_{2}^{2}} \\ {1,} & {{u_{x}^{2} + u_{y}^{2}} < k_{2}^{\prime 2}} \\ {0,} & {others} \end{matrix}\mspace{79mu} k_{i}} = {N\;{A_{i}/{\lambda\left( {{i = 0},1,2} \right)}}}} \right.$

Thus, ATF_(PC-B)(u) and PTF_(PC-B)(u) of the apodization phase contrast can be expressed as

${A\; T\;{F_{{PC} - B}(u)}} = \frac{\begin{matrix} {- {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}{P_{{PC} - B}\left( {u + u_{j}} \right)}} +} \right.}}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - B}\left( u_{j} \right)}}^{2}d^{2}u}}}$ ${P\; T\;{F_{{PC} - B}(u)}} = \frac{\begin{matrix} {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}{P_{{P\; C} - B}\left( {u + u_{j}} \right)}} -} \right.}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{PC} - B}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$

where k₀ corresponds to the radius of the objective limit, k₁ is the outer diameter of the phase contrast ring with phase delay, k₂ is the inner diameter of the phase contrast ring with phase delay. In the simulated experimental plots shown in FIG. 4(a 2), FIG. 4(a 4), FIG. 4(b 2), FIG. 4(b 4) and FIG. 4(d 1), the NA₁=0.28, NA₂=0.25, NA_(o)=0.4 are used. The phase difference ring with phase delay corresponds to the solid-filled ring in FIG. 3(d), and there is a π/2 phase delay with intensity attenuation effect, and the attenuation form is shown in FIG. 3(e). k′₁ and k′₂ is the overall outer and inner diameter of the phase contrast ring corresponding to the slash-filled ring in FIG. 3(d), there is only an intensity attenuation effect without phase delay, and the attenuation form is shown in FIG. 3(e). The algorithm can be used to adjust the width of the diagonal filling ring in FIG. 3(d), and to obtain different images of the variant phase contrast by changing the value of the numerical aperture of the objective NA₀, where the simulation results of changing the numerical aperture of the objective NA₀ can be seen in FIG. 4(d 1) to FIG. 4(d 4). The corresponding results of ATF_(PC-B)(u) and PTF_(PC-B)(u) are shown in FIG. 4(a 4) and FIG. 4(b 2). Thus, when the sample is partially illuminated by spatially coherent light, the obtained intensity spectrum of the apodization phase contrast I_(PC-B)(u) can be written as the sum of the intensity spectrum at multiple illumination angles

I _(PC-B)(u)=B _(PC-B)δ(u)+A(u)ATF_(PC-B)(u)+iΦ(u)PTF_(PC-B)(u)

This procedure results in the intensity image in the phase contrast (PC) imaging mode: the intensity image I_(PC-A)(u) obtained by the annular phase contrast method and the intensity image I_(PC-B)(u) obtained by the apodization phase contrast method, respectively.

FIG. 5 shows the experimental results of the HeLa cell, i.e., the phase contrast and differential interference contrast intensity images obtained from the quantitative phase image obtained from the TIE shown in FIG. 5(a), using the above algorithm. Among them, FIG. 5(b) displays the annular positive phase contrast intensity map; FIG. 5(c) shows the apodization positive phase contrast intensity map; FIG. 5(d) shows the DIC intensity image; FIG. 5(e) shows the annular negative phase contrast intensity image; and FIG. 5(f) apodization negative phase contrast intensity image. Comparing the phase contrast and differential interference contrast image with the quantitative phase image, it can be seen that the phase contrast and differential interference contrast image obtained by the present algorithm can achieve higher contrast and richer details; comparing the phase contrast and differential interference contrast images via the proposed method with the phase contrast and differential interference contrast images obtained by the traditional method, it can be seen that the present algorithm can achieve the same imaging effect as the traditional phase contrast and differential interference contrast images. In summary, the present invention can endow the bright-field microscope with the ability to achieve phase contrast and differential interference contrast imaging without the need for complex modification of the traditional bright-field microscope. That is, only need to use an ordinary traditional bright-field microscope without adding any complicated components. Through the PC and DIC algorithms, this method has the advantages of quantitative, high-speed, low-cost, simple structure, and less external interference. In addition, it has the same imaging performance as the phase contrast microscope and differential interference contrast microscope, which are expensive, complex-structure, and has strict environmental conditions. 

1. A microscopic imaging method of phase contrast (PC) and differential interference contrast (DIC) based on the transport of intensity equation (TIE) is characterized in that the steps are as follows: step 1, three intensity images are acquired along the optical axis, i.e., acquire three intensity images along the optical axis, which are positive defocus image I_(z-P), infocus image I₀, negative defocus image I_(z-N); the first-order axial differentiation ∂I/∂z of the intensity I can be obtained using numerical difference estimation $\frac{\partial I}{\partial z} = \frac{I_{z - P} - I_{z - N}}{dz}$ where dz is the distance between the positive defocus image I_(z-P) and the negative defocus image I_(z-N) in the axial direction z; step 2, obtain the quantitative phase image based on the TIE, that is, the quantitative phase is obtained by substituting the first-order axial differentiation ∂I/∂z obtained in step 1 into TIE; step 3, solve TIE in the differential interference contrast imaging mode, i.e., using the quantitative phase image obtained in step 2 to obtain the intensity image DIC(x,y) in the differential interference contrast imaging mode based on the DIC imaging theory; step 4, solve the phase transfer function in the phase contrast imaging mode, and solve the intensity image in the phase contrast imaging mode, i.e., use the quantitative phase image obtained in step 2 to obtain the intensity image I_(PC-A)(u) by the annular phase contrast method and the intensity image I_(PC-B)(u) using the phase contrast of apodization method based on the PC imaging theory.
 2. According to claim 1, it is characterized by the following in step 2: solve the TIE using the first-order axial differentiation ∂I/∂z obtained in step 1: ${{- k}\frac{\partial I}{\partial z}} = {\nabla{\cdot \left( {I\;{\nabla\phi}} \right)}}$ the quantitative phase ϕ can be obtained $\phi = {{- k}{{\nabla^{- 2}\nabla} \cdot \left( {I^{- 1}{\nabla{\nabla^{- 2}\frac{\partial I}{\partial z}}}} \right)}}$ where k is the wave number, k=2π/λ, λ is the wavelength, ∇ is the gradient operator, denotes dot product, ∇⁻² is the inverse Laplacian operator.
 3. According to claim 1, it is characterized by the following implementation process in step 3: (1) firstly, a beam of light passes through a Wollaston prism and divides into two coherent optical fields with Δx, Δy shear; (2) any one of them produces another phase shift of π/2, the other one remains unchanged; (3) finally, the two beams interfered and superimposed to form an intensity interferogram, that is, the final differential interference contrast image obtained; (4) in the practical case, it is not only limited to integer shear quantities, but also fractional shear quantities can be achieved, and then the solution is performed in the frequency domain: U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)] U _(Δ)(x,y)=exp(iπ/2)FFT⁻¹ {Û(u _(x) ,u _(y))exp[i2π(u _(x) Δx+u _(y) Δy)]} DIC(x,y)=|U(x,y)+U _(Δ)(x,y)|² the above method is based on Fourier domain sub-pixel panning, where Û(u_(x),u_(y))=FFT{U(x,y)} is the Fourier transform of U(x,y), (u_(x),u_(y)) is the Fourier-domain coordinate corresponding to the space-domain coordinate (x,y), whose shear Δx,Δy values are adjustable; when the shear Δx,Δy is an integer, it can be solved directly in the space domain; when the shear Δx,Δy is a fractional number, it needs to be converted to the frequency domain; meanwhile, the size of the shear Δx,Δy determines the convexity and concavity of the object in the differential interference contrast image; and the positive and negative of the shear determines the shadow direction of the object in the differential interference contrast image, and finally the intensity image DIC(x,y) in the differential interference contrast imaging mode can be obtained.
 4. According to claim 3, it is characterized by the realization process of step 3 as follows: firstly, the object optical field generated is U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)] I₀ is the focus image obtained in step 1, ϕ(x,y) is the ϕ obtained in step 2, then the light field U_(Δ)(x,y) after shearing Δx,Δy and phase shifting π/2 is U _(Δ)(x,y)=U(x+Δx,y+Δy)=√{square root over (I/(x+Δx,y+Δy))}exp{i[ϕ(x+Δx,y+Δy)+π/2]} where Δx is the amount of shear along the x-axis, Δy is the amount of shear along the y-axis, then the intensity image DIC(x,y) formed by U(x,y) and U_(Δ)(x,y) is DIC(x,y)=|U(x,y)+U _(Δ)(x,y)|²=|√{square root over (I ₀)} exp[iϕ(x,y)]+√{square root over (I/(x+Δx,y+Δy))}exp{i[ϕ(x+Δx,y+Δy)+π/2]}|² where Δx,Δy is integer N, its value is used to adjust the amount of shear and the positive and negative is used to adjust the direction of shear; DIC(x,y) is the obtained DIC intensity image; the solution for the special case of DIC imaging is as follows: {circle around (1)} for the case of uniform intensity there has I₀=1, then U(x,y)=exp[iϕ(x,y)] DIC(x,y)=|exp[iϕ(x,y)]+exp{i[ϕ(x+Δy,y+Δy)+π/2]}|² {circle around (2)} for weakly phased objects for uneven intensity, i.e., ϕ(x,y)□ 1, so U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)]=√{square root over (I ₀)}[1+iϕ(x,y)] DIC(x,y)=I ₀{1+[ϕ(x,y)−ϕ(x+Δx,y+Δy)]} {circle around (3)} for weakly phased objects for uniform intensity, thus U(x,y)=exp[iϕ(x,y)]=1+iϕ(x,y) DIC(x,y)=1+[ϕ(x,y)−ϕ(x+Δx,y+Δy)] {circle around (4)} for the case of slowly varying intensity √{square root over (I₀)}≈√{square root over (I(x+Δx,y+Δy))}, there is DIC(x,y)=|√{square root over (I(x+Δx,y+Δy))}[exp[iϕ(x,y)]+exp{i[ϕ(x+Δx,y+Δy)+π/2]}]|².
 5. According to claim 1, it is characterized in step 4: the process of obtaining the intensity image I_(PC-A)(u) obtained by the annular phase contrast method based on the phase contrast imaging theory is as follows: using annular illumination, the phase difference plate is annular, and the intensity of the object and the optical field of the object U(x,y) generated by the phase recovery are obtained by U(x,y)=√{square root over (I ₀)} exp[iϕ(x,y)] thus, I(u) = Bδ(u) + A(u)A T F(u) + i Φ(u)P T F(u) B = ∫∫L(u_(j))P(u_(j))²d²u_(j) ${A\; T\;{F(u)}} = \frac{- {\int{\int{{{L\left( u_{j} \right)}\left\lbrack {{{P^{*}\left( u_{j} \right)}{P\left( {u + u_{j}} \right)}} + {{P\left( u_{j} \right)}{P^{*}\left( {u - u_{j}} \right)}}} \right\rbrack}d^{2}u_{j}}}}}{\int{\int{{L\left( u_{j} \right)}{{P\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ ${P\; T\;{F(u)}} = \frac{\int{\int{{{L\left( u_{j} \right)}\left\lbrack {{{P^{*}\left( u_{j} \right)}{P\left( {u + u_{j}} \right)}} - {{P\left( u_{j} \right)}{P^{*}\left( {u - u_{j}} \right)}}} \right\rbrack}d^{2}u_{j}}}}{\int{\int{{L\left( u_{j} \right)}{{P\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ where I₀ is the focus image obtained in step 1, ϕ(x,y) is the phase image ϕ obtained in step 2, u denotes the frequency coordinate in Fourier space; I(u) is the intensity spectrum distribution; B indicates background items; A(u) is the amplitude spectrum distribution of the sample; ATF (u) denotes the amplitude transfer function; Φ(u) is the phase spectrum distribution of the sample; PTF(u) denotes the phase transfer function; u_(j) is the corresponding frequency shift vector caused by a single tilted illumination, j indicates different angles of illumination; L (u_(j)) denotes the intensity distribution of the illumination source in the Fourier domain; P (u) is the pupil function of the objective lens; * denotes the conjugate operation; the above expressions are general expressions for the background term, amplitude transfer function and phase transfer function under complex illumination and are applicable to all illumination modes and light pupil functions; therefore, in the case of annular phase contrast, the light source expression for the annular illumination mode L(u) is ${L(u)} = \left\{ \begin{matrix} {1,} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {0,} & {others} \end{matrix} \right.$ the pupil function under the annular phase contrast P_(PC-A)(u) is ${P_{PCA}(u)} = \left\{ {{\begin{matrix} {1,} & {k_{1}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{0}^{2}} \\ {{0.25 \times {\exp\left( {i\frac{\pi}{2}} \right)}},} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {1,} & {{u_{x}^{2} + u_{y}^{2}} < k_{2}^{2}} \\ {0,} & {others} \end{matrix}\mspace{79mu} k},\ {= {N\;{A_{i}/\lambda}\mspace{14mu}\left( {{i = 0},1,2} \right)}}} \right.$ where (u_(x),u_(y)) is the Fourier-domain coordinate corresponding to the space-domain coordinate (x,y), λ is wavelength, k₀ is the radius of the corresponding objective limit, k₁ is the outer diameter of the phase contrast ring, k2 is the inner diameter of the phase contrast ring; the integration area of the ring illumination corresponds to the phase contrast ring, inner diameter is k₂, outer diameter is k₁; then the annular phase contrast of ATF_(PC-A)(u) and PTF_(PC-A)(u) can be expressed as ${A\; T\;{F_{{PC} - A}(u)}} = \frac{- {\int{\int{{{L\left( u_{j} \right)}\left\lbrack {{{P_{{P\; C} - A}^{*}\left( u_{j} \right)}{P_{{PC} - A}\left( {u + u_{j}} \right)}} + {{P_{{P\; C} - A}\left( u_{j} \right)}{P_{{P\; C} - A}^{*}\left( {u - u_{j}} \right)}}} \right\rbrack}d^{2}u_{j}}}}}{\left( {\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - A}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}} \right)}$ ${P\; T\;{F_{{P\; C} - A}(u)}} = \frac{\begin{matrix} {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{P\; C} - A}^{*}\left( u_{j} \right)}{P_{{P\; C} - A}\left( {u + u_{j}} \right)}} -} \right.}}} \\ {\left. {{P_{{P\; C} - A}\left( u_{j} \right)}{P_{{PC} - A}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{PC} - A}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ thus, when the sample is partially illuminated by spatially coherent light, the resulting annular phase contrast intensity spectrum I_(PC-A)(u) can be written as the sum of the intensity spectra at multiple illumination angles I _(PC-A)(u)=B _(PC-A)δ(u)+A(u)ATF_(PC-A)(u)+iΦ(u)PTF_(PC-A)(u).
 6. According to claim 1, it is characterized in step 4 based on phase contrast imaging theory to obtain the intensity image I_(PC-B)(u) obtained by the apodization phase contrast method: the difference is that the former is a step filter and the latter is a trapezoidal filter, i.e., the pupil function P_(PC-B)(u) is different, expressed by the formula ${P_{{PC} - B}(u)} = \left\{ {{{\begin{matrix} {1,} & {k_{1}^{\prime 2} < {u_{x}^{2} + u_{y}^{2}} < k_{0}^{2}} \\ {{{\frac{3}{4\left( {k_{1}^{\prime} - k_{1}} \right)}u_{x}} + \left\lbrack {1 - {\frac{3}{4\left( {k_{1}^{\prime} - k_{1}} \right)}k_{1}}} \right\rbrack},} & {k_{1}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{\prime 2}} \\ {{0.25 \times {\exp\left( {i\frac{\pi}{2}} \right)}},} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {{{\frac{- 3}{4\left( {k_{2} - k_{2}^{\prime}} \right)}u_{x}} + \left\lbrack {1 + {\frac{3}{4\left( {k_{2} - k_{2}^{\prime}} \right)}k_{2}}} \right\rbrack},} & {k_{2}^{\prime 2} < {u_{X}^{2} + u_{y}^{2}} < k_{2}^{2}} \\ {1,} & {{u_{x}^{2} + u_{y}^{2}} < k_{2}^{\prime 2}} \\ {0,} & {others} \end{matrix}\mspace{79mu} k_{i}} = {N\;{A_{i}/\lambda}}}\ ,\ {i = 0},1,2} \right.$ thus, ATF_(PC-B)(u) and PTF_(PC-B)(u) of the apodization phase contrast can be expressed as ${A\; T\;{F_{{P\; C} - B}(u)}} = \frac{\begin{matrix} {- {\int{\int{L{\left( u_{j} \right)\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}{P_{{P\; C} - B}\left( {u + u_{j}} \right)}} +} \right.}}}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - B}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ ${P\; T\;{F_{{P\; C} - B}(u)}} = \frac{\begin{matrix} {\int{\int{{L\left( u_{j} \right)}\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}P_{{P\; C} - B}\left( {u + u_{j}} \right)} -} \right.}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - B}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ where k₀ corresponds to the radius of the objective limit, k₁ is the outer diameter of the phase contrast ring with phase delay, k₂ is the inner diameter of the phase contrast ring with phase delay, k′₁ and k′₂ is the overall outer and inner diameter of the phase contrast ring; when the sample is partially illuminated by spatially coherent light, the intensity spectrum of the apodization phase contrast I_(PC-B)(u) can be written as the sum of the intensity spectra at multiple illumination angles I _(PC-B)(u)=B _(PC-B)δ(u)+A(u)ATF_(PC-B)(u)+iΦ(u)PTF_(PC-B)(u) this step finally results in an intensity image in the phase contrast PC imaging mode.
 7. According to claim 5, it is characterized in step 4 based on phase contrast imaging theory to obtain the intensity image I_(PC-B)(u) obtained by the apodization phase contrast method: the difference is that the former is a step filter and the latter is a trapezoidal filter, i.e., the pupil function P_(PC-B)(u) is different, expressed by the formula ${P_{{PC} - B}(u)} = \left\{ {{{\begin{matrix} {1,} & {k_{1}^{\prime 2} < {u_{x}^{2} + u_{y}^{2}} < k_{0}^{2}} \\ {{{\frac{3}{4\left( {k_{1}^{\prime} - k_{1}} \right)}u_{x}} + \left\lbrack {1 - {\frac{3}{4\left( {k_{1}^{\prime} - k_{1}} \right)}k_{1}}} \right\rbrack},} & {k_{1}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{\prime 2}} \\ {{0.25 \times {\exp\left( {i\frac{\pi}{2}} \right)}}\ ,} & {k_{2}^{2} < {u_{x}^{2} + u_{y}^{2}} < k_{1}^{2}} \\ {{{\frac{- 3}{4\left( {k_{2} - k_{1}^{\prime}} \right)}u_{x}} + \left\lbrack {1 + {\frac{3}{4\left( {k_{2} - k_{1}^{\prime}} \right)}k_{2}}} \right\rbrack},} & {k_{2}^{\prime 2} < {u_{x}^{2} + u_{y}^{2}} < k_{2}^{2}} \\ {1,} & {{u_{x}^{2} + u_{y}^{2}} < k_{2}^{\prime 2}} \\ {0,} & {others} \end{matrix}\mspace{79mu} k_{i}} = {N\;{A_{i}/\lambda}}}\ ,\ {i = 0},1,2} \right.$ thus, ATF_(PC-B)(u) and PTF_(PC-B)(u) of the apodization phase contrast can be expressed as ${A\; T\;{F_{{P\; C} - B}(u)}} = \frac{\begin{matrix} {\int{\int{L{\left( u_{j} \right)\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}{P_{{P\; C} - B}\left( {u + u_{j}} \right)}} +} \right.}}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - B}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ ${P\; T\;{F_{{PC} - B}(u)}} = \frac{\begin{matrix} {\int{\int{L{\left( u_{j} \right)\left\lbrack {{{P_{{P\; C} - B}^{*}\left( u_{j} \right)}{P_{{P\; C} - B}\left( {u + u_{j}} \right)}} -} \right.}}}} \\ {\left. {{P_{{P\; C} - B}\left( u_{j} \right)}{P_{{P\; C} - B}^{*}\left( {u - u_{j}} \right)}} \right\rbrack d^{2}u_{j}} \end{matrix}}{\int{\int{{L\left( u_{j} \right)}{{P_{{P\; C} - B}\left( u_{j} \right)}}^{2}d^{2}u_{j}}}}$ where k₀ corresponds to the radius of the objective limit, k₁ is the outer diameter of the phase contrast ring with phase delay, k₂ is the inner diameter of the phase contrast ring with phase delay, k′₁ and k′₂ is the overall outer and inner diameter of the phase contrast ring; when the sample is partially illuminated by spatially coherent light, the intensity spectrum of the apodization phase contrast I_(PC-B)(u) can be written as the sum of the intensity spectra at multiple illumination angles I _(PC-B)(u)=B _(PC-B)δ(u)±A(u)ATF_(PC-B)(u)+iΦ(u)PTF_(PC-B)(u) this step finally results in an intensity image in the phase contrast PC imaging mode. 